Crystal Field Triplets: A New Route to Non-Fermi Liquid Physics 
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A model for crystal field triplet ground states on rare earth or actinide ions with dipolar and 
quadrupolar couplings to conduction electrons is studied for the first time with renormalization 
group methods. The quadrupolar coupling leads to a new nontrivial, non-Fermi liquid fixed point, 
which survives in an intermediate valence Anderson model. The calculated magnetic susceptibility 
displays one parameter scaling, going as T~ a (a ~ 0.4) at intermediate temperatures, reminiscent 
of the non- Fermi liquid alloy UCus-^Pda;. 
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Recent data for many "heavy fermion" Ce or U based 
compounds and alloys display diverging low-temperature 
magnetic susceptibility x{T) ~ T~ a and electronic spe- 
cific heat coefficients C e i/T = 7 (T) - T~ Q or ln(T), 

unlike those of a Fermi liquid (x(0), 7(0) constant) jjj . A 
number of theoretical scenarios have arisen to explain this 
non- Fermi lquid (NFL) physics, which broadly fall into 
two categories: (1) Localized, including models associated 
with peculiar symmetry allowed interactions between /- 
electron moments and conduction electrons [ 2 -p) or a dis- 
order induced distribution of Kondo scales [7M. Among 
the first class of models, multi-channel Kondo Hamil- 
tonians employing couplings between localized magnetic 
or orbital doublets (induced by crystalline electric field 
[CEF] splittings) and conduction electrons have been 
extensively studied (2) Extended, in which the 

NFL behavior is driven by coupling of electrons to low- 
lying modes induced by intersite /-moment interactions 
in proximity to a quantum critical point |) 10 1. ft is of 
considerable interest to sort out the applicability of these 
differing scenarios to real materials. 

Of particular interest is the alloy system UCus-aPda; 
pT| which displays NFL behavior in the range 1 < x < 
2.5 that has been described in terms of localized dis- 
ordered Kondo physics j7|,|| or Griffiths' phase theory 
[pTofl . This alloy appears to display two separate NFL 
regimes as a function of temperature: (1) An "impu- 
rity" regime independent of x for 8QK < T < 300K with 
X{T) ~ T- 1 / 3 and x"(Q,w,T) - lj- 1 ' 3 for htu > k B T 
(there appears to be negligible Q dependence to x") 
E,E3l. (2) A low temperature "lattice" regime in which 
X (T) ~ 7 (T) - T-a{x) and for which 

some evidence 

exists for intersite interaction effects The "impu- 

rity" regime data is not compatible with any plausible 
multi-channel Kondo model assuming ground state mag- 
netic or orbital doublet levels on the U ions. The picture 
is further complicated by analysis of photoemission data 
which suggests the U ions are in the mixed valent regime, 
possessing a nearly 50-50 mix of ground state weight in 
the f 3 and / 2 configurations |Q. This is problematic 
in that the Kondo effect has been conventionally studied 



only for nearly integral valent ions. 

In this paper, we present a Kondo model for uranium 
ions with a CEF triplet ground state that is allowed in 
cubic symmetry. This model features an effective spin 
1 local moment coupled via magnetic dipole and electric 
quadrupole interactions to one band of effective spin 3/2 
conduction electrons. The model displays a new stable 
NFL fixed point (FP) at low temperatures which is ro- 
bust even in the mixed valent regime of the more funda- 
mental Anderson model. We find three different power- 
law regions for the magnetic susceptibility x(T): a char- 
acteristic T~ a (a « 0.4) dependence in the intermedi- 
ate temperature region over two decades, a Curie- Weiss 
law T^ 1 at higher temperatures, and a universal power- 
law behavior T~ 2 / 3 at lower temperatures. Despite the 
differing temperature regimes, a surprising one parame- 
ter scaling emerges for x(T). The quadrupolar coupling 
is demonstrated to be relevant using multiplicative and 
numerical renormalization group (NRG) methods ||[l5) ■ 
We compare the FP properties with those of the unstable 
FP reached in the absence of the quadrupolar coupling, 
and examine the FP structure in the presence of uniaxial 
symmetry breaking fields. We argue that the intermedi- 
ate and low temperature regimes may be relevant to the 
physics of UCus^Pdj.. 

Assuming a dominantly tetravalent (5/ 2 ) U ion for the 
moment, the Hund's rule ground state has total angular 
momentum J = 4, which is split into a quadrupolar (T3) 
doublet, two magnetic triplets ^4^5) and a singlet (ri) 
under the action of the cubic CEF. The possible ground 
states, which are accessed by varying the two parameters 
of the crystal field Hamiltonian, are T3, T5, or Ti |fl6| . 
As argued elsewhere, the T3 ground state will give rise 
to the two-channel quadrupolar Kondo effect upon cou- 
pling to conduction electrons diHH- ^ ne two-channels 
arise from coupling to a local quartet (r 8 ) of conduction 
electrons which may be described as a tensor product 
of states with two magnetic labels and two orbital (r 3 ) 
labels. However, the analysis of the Tg tensors reveals 
two possible "dipole" operators: in addition to the set 
of effective magnetic doublets labeled by the orbital in- 
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dices, it can be written as a single s = 3/2 manifold. We 
find that the latter dominates the coupling to the triplet 
ground state, and this model has not been studied before. 

The f 2 triplet has both magnetic and quadrupolar mo- 
ments, and can be represented by a pseudo-spin 5=1 
which can couple both to the above mentioned Tg quar- 
tet of conduction electrons and to a local doublet (Tq 
or When the coupling to the quartet is larger, the 
corresponding Kondo Hamiltonian can be written as 
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Here the S c = 3/2 (m, m' = ±1/2, ±3/2) spin operator 
represents the four-fold T 8 states of the conduction elec- 
trons and a potential scattering term is neglected. The 
quadrupolar operators are given by {Q\ i — 1, ..,5} = 



{S y S z 
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S$, {2S 2 -Si- S 2 )/V3 }. The conduction electron with 
wave number k and pseudo-spin m has kinetic energy ek 
and is created (annihilated) by the operator c\ m (ckm)- 
In the limit of small hybridization between the conduc- 
tion band and the /-orbitals, this Kondo Hamiltonian (jl|) 
can be directly derived from an Anderson Hamiltonian 
where we restrict ourselves to the valence fluctuation be- 
tween the 5/ 1 r 7 and 5/ 2 r 4 (T 5 ) states. In this case, we 
obtain a coupling ratio Jb/ Jq — 2 and a marginally irrel- 
evant potential scattering with amplitude Jd . Our NRG 
calculations show that even including the 5/ 3 configura- 
tion and extending the Anderson model parameters to 
the rnixed-valent regime, the Kondo model ([j]) describes 
a stable FP [ jl7| ] - Therefore, at low enough temperatures, 
we can use this model to study the realistic Kondo effect 
corresponding to this FP. 

The relevance of the quadrupolar coupling Jq in 
Eq. (|l|) can be immediately seen from a multiplicative 
renormalization group procedure, valid in the weak cou- 
pling regime. After a straightforward but lengthy calcu- 
lation we derive the following RG equations: 



dj D /dx = (jrS + 12j Q )(l-5j D ) 
dj Q /dx = 6j'dJq - 36j'q - IBj^jQ 



(2) 



Here x — ln(Ep/u>) denotes the scaling variable (with Ep, 
the Fermi energy, and w, the characteristic energy scale), 
and we have introduced the dimensionless couplings jq = 
QqJq and ju = QoJd with qq, the density of states at 
the Fermi surface. The flow diagrams obtained from a 
numerical solution of Eq. (||) are shown in Fig. 1. In the 
absence of quadrupolar exchange the model scales to the 
dipolar FP "D" at (j Q = 0, j D = 1/5). This FP has been 
shown to be characterized by a critical exponent A = 1 /6 
associated with its spin sector, which can be mapped to 
the spin sector of the 10-channel Kondo problem p8| , ^9[ . 
Obviously, this dipolar FP is unstable to quadrupolar 




FIG. 1. Scaling trajectories obtained from Eq. (g). 

perturbations and for any non-zero jq it flows to a new 
FP"S" at (jd = 1/5, jq = yi760). 

At the FP " S" the ratio jq/j'd takes the value j'q/j'd = 
y^/12, and the interaction part of the Hamiltonian can 
be written in the following pscudo-S'C/ (3) invariant form: 



J C L( A c)rmn'C fc 



(3) 
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where the A"s denote the 3x3 Gell-Mann matrices sat- 
isfying the SU(3) Lie algebra [A\A J ] = 2if^ k X k , and 
can be easily expressed in terms of the spin one impu- 
rity operators. The 4x4 matrices A* are constructed 
from the conduction electron spin operators S c , and sat- 
isfy a "pseudo-S'f7(3)" Lie algebra: [A' c , A£] = 2ifi k \ k c + 
octupolar terms. While these latter terms, which arise 
from the commutators, spoil the SU(3) symmetry of the 
local triplet, they cannot couple to the impurity, and are 
hence irrelevant in a renormalization group and general 
sense. Since both the leading (second order) and next 
leading logarithmic (third order) scaling equations result 
in the FP structure of Eq. (||), we believe that this re- 
sult is universal and independent of the weak coupling 
approximation. 

Unfortunately, the pseudo-symmetry found is not 
strong enough for the usual characterization of the FP 
by boundary conformal field theory (CFT) p3| : in par- 
ticular, it is impossible to absorb the impurity spin in the 
conduction electron currents without violating the Kac- 
moody algebra of the conduction electrons. However, we 
can study the properties of the novel FP by using the 
NRG. Following Wilson (TBI, we rewrite Eq. (|lj) as 



H N+1 = A^ 2 H N + J2(fUi,rJNm + H.c), 

m 

Ho = £ £ /o : m'/o,m JT(T c ) m 'm-T 



mm' T=S,' 



(4) 



where Hq represents the effective exchange interaction 
on the impurity site, f\j m (/iVm) creates (annihilates) a 
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FIG. 2. Finite-size energy spectrum of the fixed points D 
and S of Fig. 1. The energy E is measured in units of tv/L 
(vf — 1) with L, the system size. (D) For Jq = 0, each 
energy level is labeled by axial charge q and total spin j. (S) 
For Jq 7^ 0, q is not a good quantum number since the charge 
5(7(2) symmetry is broken by the quadrupolar interaction. 

conduction electron in the logarithmic discretized band, 
and J = 2j0 O A~ 1 / 2 /(l + A -1 ) with A, the discretization 
parameter. We follow the usual procedure and iteratively 
diagonalize -ff/v+i to probe the system size on a scale of 
order fcp l \ N / 2 and temperature of order TpA~ N / 2 . 

Fig. 2 shows the finite-size energy spectrum obtained 
at FP's "D" and "S". The Jq = spectrum in the sec- 
tor with dipolar coupling only ("D") coincides with the 
exact CFT spectrum for the model of the impurity spin 
coupled to spin 1/2 electrons In the sector with a 

finite value of Jq ("S"), we show the energy spectrum we 
anticipate for a CFT of the model inferred from our NRG 
calculations (even iteration) for initial small values of Jd 
and Jq. The axial charge operator q in Fig. 2 is defined 

by q+ = E~=o(- 1 )"(/i3/2/n 1 -3/2 - 4,1/24,-1/2). and 
1z — \ E,T=o Em (fLnfnm - 5), and satisfies the SU(2) 
Lie algebra, [q z ,q±] = ±q± and [q+,q-] = 2q z Q. It 
is easy to show that the relevant quadrupolar coupling 
breaks the axial charge SU(2) symmetry down to (7(1), 
while it conserves that of the total spin j. This symmetry 
breaking, however, cannot be described by a simple phase 
shift, as is often the case: as shown in Fig. 2, some orig- 
inally degenerate SU(2) axial charge multiplets are split 
at the new FP "S" by some quadrupolar charge opera- 
tor, while others are split by a dipolar charge operator. 
These splittings are universal, apart from some trivial po- 
tential scattering. This latter is generated even for small 
Kondo couplings and it becomes more pronounced as the 
couplings are increased. 

To determine the dimension of the leading irrelevant 
operator that governs the new FP, we carried out a finite- 
size analysis of the NRG levels. Within the NRG scheme 
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FIG. 3. Top: Temperature dependent magnetic suscepti- 
bility as a function of log 10 (T). From top to bottom, the cou- 
plings are 4J D = 8 Jq = 0.5, 0.7, 1.0, 2.0, 4.0 and 8.0. Bottom: 
Rescaled magnetic susceptibility. All the data collapse to a 
universal curve. Inset: Relaxation of NRG energy levels in 
the vicinity of fixed points. The top line is for Jd = Jq = 0.1 
and the middle line is for Jd = 1.0 and Jq = 0. For compar- 
ison, we give the data for the two-channel Kondo fixed point, 
scaling as A _iY//4 . We take A = 3 here. 

the finite-size energy 1/L corresponds to ~ A~ N / 2 and 
the levels relax to their FP values as 

E NRG - £* oc K- AN '\ (5) 

where E* is the FP energy and A denotes the scaling 
dimension of the leading irrelevant operator. As shown in 
Fig. 3, for Jq ^0 almost the same energy level relaxation 
is found as for Jq = over the whole parameter space 
of Fig. 1, in agreement with an exponent A = 1/6. This 
implies that \{T) and 7 (T) behave like T 2 ^ 1 = T~ 2 / 3 
at low temperatures. Since the S = 1 impurity spin also 
incorporates orbital degrees of freedom, similar behavior 
is expected in the orbital (stress) susceptibility, or the 
stress-induced magnetization. 

To determine the impurity susceptibility x(T) we cal- 
culated the temperature dependent magnetization in- 
duced by a small local field at the impurity site. The 
resulting curves are plotted in Fig. 3. Each curve has 
an interesting region where it behaves like T~ a over ap- 
proximately two decades, and 1/3 < a w 0.4 < 2/3 
slightly depending on the magnitude of the Kondo cou- 
plings. After adjusting the overall scale of x(^)j as shown 
in Fig. 3, we can place all the x(^) data on a single uni- 
versal curve using a single temperature scale T* (though 
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strong coupling induces some deviations from scaling at 
higher temperatures-see the discussion below) . This one 
parameter scaling in temperature strongly suggests that 
the intermediate temperature regime behavior reflects 
the new low temperature FP rather than some unstable 
FP. For small to intermediate Kondo couplings x behaves 
according to the Curie- Weiss law T~ a (a = 1) at large 
temperatures. For smaller temperatures an intermediate 
region appears where a ~ 0.4. As the temperature de- 
creases further, x(T) turns up and behaves as ~ y- 2 / 3 
in the vicinity of the novel FP. For larger couplings the 
Curie-Weiss part is absent, and x starts as ~ T -1 / 3 at 
high temperatures and then the exponent a gradually 
goes up to 2/3 at low T. The appearance of the in- 
termediate region with x{ T ) ~ T ~ a (1/3 < a < 2/3) 
is specific to non-zero quadrupolar coupling Jq. When 
Jq = 0, the dipolar coupling Jd gives only monotonic 
behavior T^ 1 — > T~ 2 / 3 with decreasing temperatures. 
This is also a clear difference between the Kondo effect 
for local triplet and doublet states. 

Finally, we discuss the stability of the novel FP against 
a uniaxial (tetragonal) lattice distortion. The distortion 
lifts the triplet degeneracy, giving either a singlet or dou- 
blet ground state split by a value /iq, and destabilizes 
the novel triplet FP. Instead, we have two possible sta- 
ble FP's: when the singlet lies lowest, a Fermi Liquid 
FP arises, and a ground doublet experiences a NFL FP, 
associated with the two-channel quadrupolar Kondo ef- 
fect for tetragonal symmetry H . According to our NRG 
calculation, the crossover temperature below which the 
system flows away from the cubic FP to one of the above 
two varies as h,Q for a small distortion |l7| ]. This power 
implies that the operator with j = 2 corresponding to 
lattice distortions has the same scaling dimension 1/2 
at "S" as at the unstable FP "D" of Fig. 1. The two 
tetragonal FP's are separated by a boundary line in cou- 
pling space on which the novel cubic FP resides. The 
features associated with this cubic symmetric FP are ex- 
pected to appear in some U based compounds with uni- 
axial anisotropy, too. 

In conclusion, we investigated a new Kondo Hamil- 
tonian describing the dynamics of a local triplet. The 
model is not restricted to uranium ions alone, and could 
possibly describe praseodymium, thulium, and terbium 
based materials as well. The quadrupolar exchange inter- 
action drives the model to a new fixed point, character- 
ized by a universal charge SU(2) symmetry breaking, a 
leading irrelevant operator with dimension A = 1/6, and 
a pseudo-S*[/ (3) symmetry. The universal magnetic sus- 
ceptibility has an intermediate temperature range, where 
it scales as x ~ T ~ a wit h 1/3 < a « 0.4 < 2/3, with 
strong coupling producing a reduction of the power law 
towards a = 1/3. Since the extreme mixed valence of 
UCu5_ x Pd x suggested by photoemission is compatible 
with a strong coupling limit of the Kondo model, we sug- 
gest that the intermediate temperature range susceptibil- 



ity may correspond to the "impurity" range identified for 
this material. The surprising increase of the power law 
at lower temperatures will give the concentrated system 
a greater tendency towards intersite interaction effects, 
qualitatively compatible with the suggested interaction 
driven low temperature physics. To test the idea further, 
we strongly urge an experimental study of this system 
with uranium diluted away by thorium. The model may 
prove relevant to the alloy Yi-^U^Pds, too, for which 
recent neutron data suggests nearly degenerate T3 and 
r§ states on the uranium ions |22| . 
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